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ABSTRACT 

We develop a perturbative approach to study the excitation of disturbances in the 
extended atomic hydrogen (HI) discs of galaxies produced by passing dark matter 
subhaloes. The shallow gravitational potential of the dark matter subhaloes (compared 
to the primary halo) allows us to use the epicyclic approximation, the equations of 
which we solve by modal analysis, i.e., assuming a disc is composed of N radial rings 
with M modes. We show that properties of dark matter subhaloes can be inferred from 
the profile and amplitude of the modal energy of the disc. Namely, we find that the 
overall amplitude of the response gives the mass of the dark sub-halo. Motivated by 
this modal analysis, we then show that the density response shows similar features. 
Finally, we show that our results agree with those from full hydrodynamic simulation. 
We find a simple scaling relation between the satellite mass and Fourier amplitudes 
of the resultant surface density of the gas disc where the effective Fourier amplitude 

1 /2 

(essentially a sum over the low order modes) scales as nis , where is the satellite 
mass. The utility of this relation is that it can be readily applied to an observed HI map 
to deduce the satellite mass without recourse to full numerical simulation. This will 
greatly aid us in analyzing large samples of spiral galaxies to constrain the population 
of dwarf satellites in the Local Volume. 

Key words: galaxies: interactions - dark matter - galaxies: general 



1 INTRODUCTION 

Dark matter halos grow by a sequence of mergers and ac- 
cretion events in cold dark matter (CDM) cosmology. Due 
to the absence of scale in this heirarchical structure forma- 
tion scenario, dark matter haloes are expected to be merely 
scaled-up or down versions o f each other with the total mass 
being the scaling parameter (|Navarro et al.lll997l '). Observa- 
tion of galajcies at different masses show that this is clearly 
not the case. At the high end, clusters have thousands of 
;alaxies, i.e., subhaloes, in li ne with theoretical predictions 
Nataraian fc SpringeJ [20041 ). However, galactic size halos 
have a few bright satellites, e.g., our Milky Way (MW), 
in constrast with th eoretical models of the MW halo (e.g. 
iDiemand et al.|[2008l ). and dwarfs have no bright satellites. 
This qualitative difference at the faint end has come to 
be k nown as t he substr ucture or "missing sa t ellites" prob- 
lem llKlvpin et al. 19£9l: iKravtsov et al.|[20oi : iMadau et all 
l2008l : lKravtsovll201C 



* E-mail: pchang@cita.utoronto.ca (PC); schakral@fau.edu (SC) 
^ To be precise, the "missing satellites" problem refers to the 
lack of observed satellites in the MW halo compared to what is 
expected theoretically. 



iKravtsovl (|2010t l recently review the possible solutions 
to the "missing satellites" problem. To summarized, solv- 
ing the "missing satellites" problem demands that a scale 
be introduced either in structure formation of dark mat- 
ter haloes or in the physics of galaxy formation to suppress 
the number of observed systems at small scales. Namely, 
the power spectrum of fluctuations must be suppressed at 
small scales via warm dark m atter ( Zcntncr & Bullock 2003: 
iHooper etTaribOOTl : |Primackl l2009') or s tar formation is sup- 



pressed in dwarfs due to reioniza tion (Bul lock et al 



lOkamoto fc Frenkll2009l :[ Busha et al. 2010; Iliev et al 
or feedback from the first stars (Iwise fc Abeill2008l ). 



I200C : 



20101 ) 



Recent discoveries of dwarf galaxies in the MW and 
M31 have rel ieve some of this tension (for a review see 
IWillm an"2010'). These dwarfs are the most dark matter dom- 
inated ( Strigari et al. 200§) and metal poor systems known 
(|Kirbv et al.ll20 08). While their unusual properties point to 
the difficulty of forming stars in these systems, the num- 
ber density of these dwarfs may alleviate the "missing satel- 
lite" problem_once completeness corrections are accounted 
for (|Tollerud et al.|[2008 ). 

The evidence from the MW observations and theoretical 
predictions point to the existence of many nearly da r k sub - 
haloes hosting few stars. For instance, iBullock et al] (|2010l ) 
recently pointed out that these stealth galaxies would be 



2 Chang & Chakrabarti 



missed in current survey, having surface brightnesses that 
are too diffuse. This motivates finding alternative means to 

constrain this dark subhalo populat ion. 

In a recent series of papers (|Chakrabarti fc Blit j|2009l . 

I2OI0I . hereafter CB09 and CBIO, respectively), CB09 and 
CBIO made the assumption that the observed distur- 
bances in the oute r atomic hydrogen gas disc of the MW 
l|Levine et al.i '20061 are due to the tidal interaction with a 
dark matter dominated dwarf galaxy. Proceeding under this 
hypothesis, CB09 and CBIO numerically simulated a suite 
of encounters with dwarf galaxies interacting with the MW, 
and developed an approach to infer the mass and current 
location of satellites from analysis of their tidal imprints on 
outer gas discs. The Fourier amplitudes (CB09) and phase 
(CBIO) of the suite of simulations were compared to the 
observed amplitudes and phases to determine the pericen- 
ter approach, current position (in radius and azimuth), and 
mass of the subhalo that drove the observed disturbances. In 
doing so, CB09 and CBIO demonstrated how to characterize 
a dark subhalo from looking at it's effects on the primary 
galaxy's extended HI disc. Following this work, we demon- 
strated the validity of this method by performing the analy- 
sis with respect to galax ies with known optical companions 
ijChakrabarti et al.ll201ll , hereafter CBCB). We showed that 
the method accurately recovers the mass of the satellite and 
its position purely from analysis of observed disturbances in 
the HI disc, without requiring any knowledge of the optical 
light from the satellites. 

In this paper, we continue this work and develop a sim- 
plified test particle approach to studying the generation of 
these disturbances in extended HI discs by dark matter sub- 
haloes. Using this approach, we develop scaling relations be- 
tween the satellite mass and Fourier amplitudes of the re- 
sultant surface density of the HI disc. These relations can 
be utilized by observers who wish to determine the satellite 
mass directly from the observed HI map, without having to 
take recourse to full numerical simulations. Our simplified 
approach is motivated by the fact that DM subhaloes have a 
much lower velocity dispersion, i.e., gravitational potential, 
than the host halo in which they sit. Thus, we calculate the 
motions of the gas (or stars) in the epicyclic approximation 
and show that this approximation matches the results of 
both test body calculations and full numerical simulations. 

This work is complementary to the work of CB09, 
CBIO, and CBCB. Rather than compare the results of nu- 
merical simulations to observations of specific extended HI 
discs, this work focuses on the physics of the subhalo exci- 
tation and exploration of its parameter space. In the future, 
we will apply these scaling relationships to a large sample of 
local spiral galaxies to determine the statistical viablity of 
the Tidal Analysis method, i.e., the incidence of false posi- 
tives. Moreover, we will attempt to constrain the population 
of dwarf galaxies in the Local Volume. 

This paper is organized as follows. In iJH we take the 
equations of motion for a test particle and develop its modal 
equivalent in the epicyclic approximation. We then show 
how the modal viewpoint can be related to the test particle 
viewpoint. This modal analysis is shown to reproduce test 
particle calculations accurately. In addition, the invariance 
of the modal energy response (in the absence of dissipation) 
provides us an unique vantage point to study the response 
of the disc. In ^ we study the modal energy response and 



the corresponding density response. We also compute global 
quantities from these responses as a measure of the global 
properties of the disturbances raised by the passing DM sub- 
halo in ij4] We demonstrate a scaling relation between the 

disc response and the mass of the satellite where the ampli- 

1/2 

tude of the modes scale like Ms . To show that our analysis 
accurately models the relevant physics, we then compare our 
results to those computed from large scale simulations in ^ 
The good agreement between the SPH simulations and our 
results is encouraging. In addition, the scaling which we de- 
rived from our simplified calculation continues to hold. Fi- 
nally, we close in ^Hlwith a discussion of how these methods 
can be used to constrain the properties of DM subhaloes in 
galaxy haloes. 



2 BASIC EQUATIONS 

We begin with the equations of motion of a test particle in 
2-D {r,6 coordinates): 
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where is the radial velocity, J — 4>r'^ is the angular mo- 
mentum, <jf>(t) is the 6 position of the particle, and "I> is the 
gravitational potential. We assume a background state of 
circular orbits, i.e., 11^,0 = 0, about an axissymmetric poten- 
tial, $0. Expanding equations ((T| - Q to linear order, e.g., 
r ^ r + Sr, we find 
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We assume that the perturbation to the potential, 5$, arises 
from external perturbers ignoring the self-gravity of the par- 
ticles. This simplification is safe for typical HI disc parame- 
ters. 

Assuming solutions of the form exp{im{9 — fl{r)t)), 
where fl{r) — (j>{r), e.g. 

M 

Sr (X Srm.{r,t) exp{im{9 — Q{r)t)), 

m— 1 

reduces equations (O - ([HI) to M separate wave equations, 
sourced by the Fourier components of the perturbing poten- 
tial: 
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(9) 
(10) 



- — V - -Q^ (11) 
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imS^ exp {imQ{r)t) . 



(12) 



dSJ _ d£ 
dt ^ dr 

Solving these perturbed equations for a finite M (say M — 
10) allows us to estimate the hnear response of a circular 
ring of of orbiting particles to external disturbances, i.e., 
passing DM subhaloes. Thus, any individual particle's per- 
turbed orbit, whose guiding center is at radius r, can be 
found by summing over this response. Namely, we recon- 
struct the particle positions from these modes using: 



ri{t) = ri,o -I- Sri(t), 
(j}i{t) = (j}i,o + S(jii{t), 



(13) 
(14) 



where Vi^ and (j>i,o are the initial radius and azimult of the 
particle, i, and 

M 

Sri{t) = ^ Sr,n{rifi,t) exp{im{(j)ifi -^l{rifi)t)), (15) 

m=l 
M 

5(t>i{t) = ^ 5<^m(n,o, t) exp{im{(j)ifi - n(ri,o)t)), (16) 



This greatly simplifies this orbit into a sum of M simple 
harmonic oscillators, whose natural frequency only depends 
on their radial position. 

To illustrate the power of this formalism, we compare 
a direct numerical calculation of a coUisionless particle disc 
with that of the modal calculation above and show the re- 
sults of the direct test particle calculation and our modal 
reconstruction in Figure [T] Here, the broad features of the 
response is captured between the mode calculation and the 
N-body calculation. It is evident that the modal calculation 
captures the details of the perturbation accurately. 

Our formalism makes the analogy between the differ- 
ent modes of a radial ring and a collection of uncoupled 
simple harmonic oscillators, whose oscillator frequency just 
depends on their radius. Hence it is natural to think of the 
modal energy, 



E„ 



(17) 



where k = ^ R{dQ? /dr + AO? is the epicyclic frequency and 
20,5, J m 



5r' = 5r„ 



(18) 



is the perturbed radial mode m that accounts for the per- 
turbation to the guiding center. Looking at the oscillators 
energy has the major advantage that it remains constant in 
the absence of dissipation or excitation, whereas the density 
response of a disc varies with time. This is advantageous 
in elucidating the physics. In the next section, we show the 
benefit of this mode of analysis and relate its properties to 
the density response. 



3 DISC RESPONSE 

For the purposes of this work, we make a series of simpli- 
fying assumptions. First, we assume a static potential for 
both the primary galaxy and the subhalo, i.e., we do not 
use a "live halo". Second, we assume a constant density 
HI disc motiv a ted b y observations of the Galactic HI disc 
IWong fc Blit3 ll200dl and the extended HI discs of external 
galaxies iBigiel et al. (|2010l ). There also exists an exponential 



gas disc, but we ignore this for the sake of simplicity0 These 
assumptions allow us to develop a simplified approach to the 
modeling of disturbances in extended HI discs. We demon- 
strate the accuracy of the results of the simplified approach 
to that of a full simulation later in iJS] 



3.1 Modal Energy Response 

As we mentioned earlier in 321 the modal equations reduce 
to a number of simple harmonic oscillators. A passing per- 
turber excites these simple harmonic oscillators, each to a 
different amplitude and phase. In the absence of dissipation, 
the modal amplitude or energy remains constant. We illus- 
trate this point in Figure [2] where we show six snapshots of 
the orbit of a coplanar subhalo with rpcri = 20 kpc. Initially, 
the energy is zero, but rises dramatically after the disc suf- 
fers an encounter, before asymptoting to a constant. This 
property of the modal energy response makes them illumi- 
nating for studying the response of the disc. We should point 
out that while the energy is already substantial when the 
perturber is at pericenter, the density response of the disc 
does not happen until much later. This is because, the en- 
ergy imparted by the perturber at pericenter serves mainly 
to increase the perturbed velocity. However, this perturbed 
velocity only translates into perturbed density after ~ dy- 
namical time. 

How does the modal (energy) profile depend on the or- 
bital parameter? To answer this question, let us first fix 
the mass and the orbit to be that of a parabolic {E = 0) 
coplanar orbit so that the modal profile only depends on the 
pericenter distance of the perturber. In Figure [S] we plot 
modal energies as a function of radius for a 1:100 encounter 
for r-peri = 10 kpc (left plot) and 20 kpc (right plot). Note 
the dramatic increase in modal energy for the oscillator at 
r = rperi. Indeed this increase in modal energy is seen for 
1:10 and 1:1000 coplanar encounters as shown in Figure U) 
We also note that the amplitude of the modal profile de- 
pends only on the mass and is only weakly dependent on 
Tpcri at least at large distances. 



We now fix rn 



20 kpc to study the effects of in- 



clination in Figure (5] In the coplanar case studied above, 
all the modes showed a similar feature at rpcri- Here, we 
instead see that only the m = mode consistently shows 
this same feature as the coplanar case. Specifically, for all 
inclinations, ranging from +T7/2 (prograde coplanar) to 
vertical to —n/2 (retrograde coplanar), the m = modal 
energy increases dramatically as r > rpcri. 

The effect of inclination on the m 7^ modal energies 
is more complicated. For the coplanar prograde interaction, 
we also see a large increase in the modal energy at r = 
rpcri for all m ^ modes. However, this is now true for 
the other interactions. For instance, the retrograde coplanar 
interaction {~n/2) shows a large amplitude for the m — 1 
mode inside of rpcri and it decreases outward. Nevertheless, 
it is clear the modal energies does show a sudden change at 

rpcri . 

In these examples above, we have studied the effect of 
varying the inclination for purely prograde encounter, i.e.. 



2 CB09, CBIO, and CBCB use an exponential gas disc with an 
added flat HI disc. 
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Figure 1. Comparison between a direct integration (left) and the mode calculation (right) of a coUisionless particle disc excited by a 
passing DM subhalo. In this calculation, we use 50000 particles for the test-particle calculation. For the modal calculation, we set M = 64 
for 200 rings. From this plot it is clear that modal and N-body calculation agree. 



the case where the disc crossing radius is equal to the peri- 
center distance. Though we do not present them here, we 
have also studied the case of fixing r-pcri and the inchna- 
tion and varying the angle at which the perturber crosses 
the disc, i.e., 7^ for the disc crossing point on the y- 
axis. In these examples, the modal energies also show a 
change at rpcri, which suggest that looking for changes in 
the modal energy in the disc is a good way to deduce rpcri- 
We briefly comment on this in the next subsection, but re- 
serve a more extensive discussion on the operational details 
to future work. 



3.2 Density Response 

The modal energies provide a good proxy for the response of 
the disc, but our analysis of modal energies here is new. We 
now link this modal energy response to the density response, 
which is the observed manifestation of these perturbations. 
To calculate the density response from the modal energies, 
we reconstruct the position of the particles using equations 
(|13|) - (|14|) . From this, we obtain the density of particles as 
a function of r, (p, and t, from which the Fourier amplitudes 
can be calculated. 

Our analysis of modal energies show that they are lo- 
calized to certain regions of the disc. This suggests that the 
density perturbations also follow a similar pattern. We show 
this in Figure [6] where we show density perturbation from 
modal theory for the same six snapshots of the orbit of a 
coplanar subhalo with rpcri = 20 kpc as in Figure [S] Note 
that we have only not included the m — mode here in 
constrast to the modal energy comparison. While, we can 
define such an axisymmetric perturbation based on the ini- 
tial (unperturbed) configuration, such a construction would 
not be well defined if we compare this to observations or 



simulations (see As we anticipate such a comparison, 
we have ignored the perturbed m — mode in the density 
response. In addition, we have normalized the m = 1 — 4 
modes relative to the complete m = mode, which is dom- 
inated by the unperturbed component to further facilitate 
comparisons to observations and simulations. 

We also show the same six snapshots using a full test- 
particle calculation in Figure[7l Note the similarity of Figure 
[6]and[7l which increases our confidence in our method of so- 
lution. Comparing Figure [2] and [H] side by side, it is clear 
that the modal energies are already significant as the per- 
turber crossed rpcri, while the density perturbations occur 
later. As discussed above, this is because the modal energies 
is initially concentrated in the perturbed velocities and the 
densities response only arises after ~ dynamical time. 

The abrupt increase in the modal energies at r = rpcri, 
which was previously identified by looking at modal energies, 
is also seen in the density response. The density response like 



the modal energies remains small in the regions r ; 



ri, but 



increase rapidly for r > rpcri- What is also interesting is that 
the identification of rpcri from the density response in Figure 
[6] and [7| seems to be clearer than the modal energy response 
plotted in Figure[2l In addition, the density response appears 
to be much more localized than the modal energy response 
(Fig.[2|. Namely, it is strong between rpcri and rpcri + 10 kpc, 
whereas the modal energy response is strong all the way out 
to the end of the disc. 

To help illiustrate the effect of this density response, we 
define a synthetic measure of the power of these modes as: 



atot(r) 



\ ^ 5 ' 



(19) 



and plot this in Figures |6] and [T] 
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Figure 2. Modal energies of tlie disc in response to the disturbances by a coplanar 1:100 perturber witli an Tpcri = 20 kpc. In tlie upper 
left hand corner of each plot is the radial distance of the perturber. Initially when the perturber is far away (a), only the other disc 
responds to the perturber. However, as the perturber gets closer, more and more of the disc is perturbed. In (c), the perturber is at rperi 
and we note a kink in the m = modal energy at r = r-pej;. This kink on the modal energy shows up in the subsequent plots (d), (e), 
and (f) in the other modes as well. This kink is the result on an increase modal energy response at rperi and serves as an indicator of 
'■peri of perturbing subhaloes. We also note the advantage of studying the modal energy as modal energy is fixed after the interaction 
with the perturber in the absence of dissipation. 



We now use this synthetic measure (eq.[T9]) to show localized response, especially for the m = mode allows us 



that this localized response around rperi continues to hold 
as we vary the orbit with different inclinations and different 
angles. In Figure [8] we show the effects of different inclina- 
tion on the Fourier power of the resulting density structure, 
while holding angle fixed. Comparing Figure [S] with O we 
find that again that the rapid increase in the m — power 
is a good discriminant for rperi- The m 7^ modes also show 
good discriminating power though not as the m = mode. 
However, what is really encouraging is that, while the m = 1 
mode varies signficantly with inclination in Figure O its ef- 
fect on the density response is much milder. From looking 
at these plots, it is clear that the m — modes has the most 
discriminating power, but atot appears to be adequate for 
determining rperi- 



4 SCALING RELATIONS: INFERRING 
SATELLITE MASS FROM FOURIER 
AMPLITUDES 

Our discussion in § [3] shows while the overall dynamical re- 
sponse of the disc depends on the orbital parameters, i.e-, 
rperi, inclination and angle, the density response of the disc 
is generally localized around rperi and rperi + 10 kpc. This 



to constrain rperi of the orbit. We now discuss the effect of 
the subhalo's mass on the disc response with an aim toward 
developing scaling relations between the response and the 
perturber's mass. 

We begin by defining an appropriate global measure of 
the response of the disc in the much the same spirit as our 
earlier definition of atot (eq.[T9]). As we have seen already 
that the disc response is localized, we define the effective 
amplitude of the disc as 



efl = (Ar) 



i{r)\dr 



(20) 



where Ar = 10 kpc, motivated by the results of the previous 
section. We then define the total effective amplitude as 



1 



it 



am.cff 



(21) 



A major advantage of defining such global quantities is that 
like the modal energies, their overall amplitudes settles into 
a relatively constant range after the initial excitation, in the 
absence of dissipation. We show this in Figure [51 for DM 
perturbers with mass ratios of 1:10, 1:100, and 1:1000. Note 
that at.eff settles into a constant range that depends on the 
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Figure 3. Modal energies of the disc for a 1:100 perturber with rp^ri = 20 kpc (a) and 30 kpc (b). In both case the perturber sits at 
roughly rp 80 kpc after its interaction. As the modal energy is fixed at this point in the absence of dissipation, we can look in detail 
at the structure of the interaction. Note that in both cases, a kink on the modal energy appears at r = rp^ri. This strengthens the point 
made in Figure [2] that these kinks serves as a good proxy for Tperi. 




mass ratio ~ 1 Gyr after their interaction, which is roughly 
one rotation period of the outer disc. 

Motivated by this result, we average at, off after ~ one 
rotation period for different inclinations and angles in Figure 
1101 The error bars define the variation in at,eff for different 
angles. The bands represent the variation of Ot,eff for each 
mass ratio as we vary the inclination angle from 90 to - 



90. The variation in the inclinations is « 2. However, the 
average amplitude over which at, off varies depends strongly 
on the mass ratio of the satellite. We show this mass ratio 
dependence more clearly in Figure [11] where we collate the 
results of Figure [TO] and plot the overall disc response as a 
function of tUs- The error bars represent the variation over 
inclination and angle (same as the banded structures in Fig 
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Figure 5. Modal energy profiles of a 1:100 perturber with rpcri = 20 kpc for different inclinations: i = 90 (a), 60 (b), 30 (c), (d), -30 
(e), -60 (f), -90 (g) degrees. The kink in the modal energy profile in the coplanar case (a) and Figure[2]is only persistent for the m = 
mode for all inclinations. However, the modal energy profile near r = rpcri shows features in all cases like peaks, e.g., m = I case for -30 
(e), -60 (f), and -90 (g) degrees. 
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^- -1.5 
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nop . In addition we plot the following fitting formula (dotted 
line) 



-2.5 



0.5 1 1.5 

t [Gyrs] 

Figure 9. Plot of the synthetic Fourier amplitude of the den- 
sity response as a function of time. Note that after the initial 
interaction, the amplitude stay relatively constant in the absence 
of gaseous dissipation. Whereas this is expected for the modal 
energies, it is surprising that this is preserved in the density re- 
sponse as well. It points out that the synthetic amplitude is a 
good proxy for measuring the global properties of the response in 
a robust manner. 



dt.off 



0.5 



(22) 



Equation (|22p provides a reasonable fit to the numerical re- 
sults as shown by the dashed line in Figure 1111 It suggests 
that the energy of the modes, E (x oc Ms ■ This is different 
from the expectations of the imp ulse approxi mation which 
instead suggest that the E cx (|Whitdl200Gl '). The fact our 
results differ from the impulse approximation means that 
the impulse approximation is not a good approximation for 
the interaction of subhaloes with extended HI discs, a point 
already discussed by CB09, CBIO, and CBCB. 



5 COMPARISON TO SPH SIMULATIONS 

To confirm our analysis, we compare our results from the 
modal analysis above to full-scale numerical galactic sim- 
ulations. We use the sm ooth particle hydrodynamic code, 
Gadget-2 (|Springelll2005l ) using the methodology presented 
in CB09. GADGET-2 uses an N-body method to follow the 
evolution of the coUionsionless components, and SPH to fol- 
low the gaseous component. The simulations of disc galax- 
ies tidally interacting with dark matter mini-halos reported 
here (unless otherwise noted) have gravitational softening 
lengths of 100 pc for the gas and stars, and 200 pc for the 
halo. The number of gas, stellar, and halo particles in the 
primary galaxy are 4 x 10^ , 4 x 10^ and 1.2 x 10^ respectively 
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Figure 6. Fourier amplitudes of the density response of a disc suffering a 1:100 perturber encounter with rpcri = 20 kpc. The perturber 
has a coplanar orbit (i = 90). This corresponds to the modal energy calculation of Figure[2] where we showed that the encounter generates 
a kink in the modal energy at r = rperi- Here we see that the result of this modal energy kink is reflected in the density response of the 
disc in (d), (e), and (f). Namely, the density response is peaked at r fperi- The response of m = 1 — 4 is compared to the m = mode. 
For the m = line shown, we have only included the perturbation and compared it to the initial uniform background. 



for our fiducial case. We refer the interested reader to CB09 
for additional details on the simulation methodology. 

We focus here on a few representative cases that span a 
reasonable range in parameter space. We list in Table 1 the 
parameters for the simulations that we use for this compar- 
ison. They vary in mass ratio between 1:10 and 1:300 and 
have coplanar and polar inclinations. 

As an example, Figure [12] shows the density response 
for a few time snapshots of the interaction for the 1:100 
coplanar example. The general response of the disc in the 
full simulation is visually similar our modal calculations of 
Figure[T]especially that oit = 0.650 Gyrs. This is somewhat 
unsurprising as the dominant physics in both calculation is 
the gravity of the primary halo and the subhalo. However, 
it does highlight that much of the physics that has been 
ignored in the simplified calculation, i.e., a live dark mat- 
ter halo, tidal stripping of the subhalo, gas dissipation, star 
formation, and self-gravity, does lead to some diff'erence in 
the appearance of the disturbances in the HI disc. A main 
physical effect that has not been taken into account in these 
calculations is gas dissipation. Neglecting this effect allows 
for simplicity of calculation as the modal energy (and the 
Fourier amplitudes once generated) are constant in the ab- 
sence of dissipation. In the presence of gas dissipation, the 
Fourier amplitudes decrease as a function of time as distur- 
bances in the gas disc damp out on the order of a dynami- 



cal time. We do not model this effect here, which would be 
needed to identify the current location of a satellite from 
observed disturbances in the HI disc as done earlier in the 
SPH calculations by CBCB. 

Motivated by the visual agreement of the SPH and 
modal calculation, we focus on comparing at,cff for the SPH 
simulation and the modal calculation. Here the agreement 
is striking. In the final column of Table 1, we list the at, oft 
for the various SPH simulations that we have performed. We 
also at, off for the SPH simulations as x's in Figure [TOl and [TT] 
The SPH simulation and the modal calculation show strik- 
ing agreement especially for the 1:100 case. We have also 
computed the 1:300 case and it too shows excellent agree- 
ment between the modal calculation and the SPH calcula- 
tion. The 1:10 case show some mild disagreement between 
the SPH simulation and the simplified modal calculation. 
This is not unexpected as linear theory should start to break 
down for sufficiently massive satellites (mass ratios > 1:10). 
The overlaid points from the SPH simulation in Figure 1111 

1/2 

gives additional support to the at, eft cx: rris scaling that we 
initially found from the modal calculation. 



6 CONCLUSION AND DISCUSSION 

Subhalo interactions with the HI discs of galaxies should 
leave measurable effects on the gas long after the subhalo 
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Figure 8. Snapshots of 1:100 perturber rperi = 20 kpc for different inclinations. For each inclination, we have taken the same snapshot, 
which corresponds to (f ) of Figurc|6] The density response peaks around r = rperi in each of these snapshots, in Hne with the expectations 
from our energy analysis. 
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Figure 12. The time evolution of the gaseous density images of a Milky Way like galaxy interacting with a 1:100 dark sub-halo with a 
pericentric approach distance of 33 kpc. The units of time, marked in the figure, are in Gyr. The box extends from -60 kpc to 60 kpc. 




100 



Figure 10. Late time amplitude as a function of the inclination 
for 1:10 (top set of points), 1:100 (middle set of points), and 
1:1000 (bottom set of points). The error bar represent the l-tr 
spread in the synthetic amplitude as the angle is varied and the 
banded regions represents the l-cr spread in the range of synthetic 
amplitude as inclination is varied. 




Figure 11. The central and the ranges of at eff ^ a- function of 
mass ratio, varying from 1:1000 (left point) to 1:10 (right points). 
The x's show the results from the full SPH simulations described 
in ^ The dotted line shows a relation that scales like rris''^. 



has passed by as demonstrated by CB09, CBIO, and CBCB. 
In this paper, we look carefully at the dynamics of these 
interactions in the context of a simplified epicyclic approxi- 
mation. 

• By performing a modal analysis on the equations of 



motion, we show that we can reproduce the orbital dynamics 
and distribution of particles of a test-particle calculation 
for identical conditions. While this may not seem surprising 
given the mass range of the perturbers that we are interested 
in (< 0.1), we show that performing such an analysis, we can 
understand crucial aspects of the nature of such interations. 
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Table 1. Parameters of SPH simulations and ae,tot 



mass ratio 


rpcri [kpc] 


orientation 


Cc.tot 


1:10 


30 


coplanar 


0.22 






polar 


0.11 


1:100 


32 


coplanar 


0.061 






polar 


0.055 


1:300 


34 


coplanar 


0.039 






polar 


0.036 



In particular, we show that the modal energy shows sharp 
variations at ~ rperi- 

• By developing a sense of what to look for from the 
modal analysis, we then showed that such variations also 
show up in the Fourier modes of the perturbed density. 
Namely, we find sudden increases in the Fourier modes of 
the density response at rpcri- Such sharp variations at 
large radii are a distinct signature of tidal interaction and 
would argue against a secular origin. 

• We also showed that this simplified analysis agrees 
with test particle calculations as well as more extensive SPH 
simulations. This demonstrates that the physics ignored in 
the simplified approach (a live dark matter halo, tidal strip- 
ping of the subhalo, gas dissipation effects, self gravity and 
star formation, and dynamical friction) do not significantly 
alter the production of the disturbances in the HI discs. 

• One of our main results here that can be readily em- 
ployed by observers is that the effective Fourier amplitdues 
of the resultant surface density scales as the square root of 

1/2 

the mass of the satellite, i.e., at, off oc ms . This suggests 
that the energy of the modes is cx nis. This scaling rela- 
tion can be used to analyze observed HI maps and detemine 
the mass of the perturbing satellite without having to take 
recourse to full numerical simulation. 

In future work, we will apply these scaling relations to 
analyze results from cosmological simulations to determine 
how the dark sub-halos will impact the disc. To do this, 
we will generalize our method to incorporate multiple per- 
turbers. The effect of multiple perturbers is not yet known, 
but it may be the case that multiple perturbers do not sig- 
nificantly alter our results. Namely, if the timescale between 
encounters is significantly longer than a dynamical time, the 
disturbances in the gas disc from the first encounter will be 
damped out. This short — term memory of the gaseous com- 
ponent allows us to more cleanly disentangle the effect of the 
last perturber from the ones that impacted the disc previ- 
ously. Simulations predict impacts with '-^ 1 : 100 mass ratio 
pertu rbers occur only once every ~ Gyr (|Kazantzidis et al] 
I2OO8I ). which bodes well for this possibility. 
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